CO 


LANGLEY RESEARCH CENTER 


76 00507 2 


95 




NASA Contractor Report 165847 


NASA-CR-165847 

0- OG I \ 2-^ ^ 

A VORTEX-FILAMENT AND CORE MODEL FOR WINGS ~ ^ ^ 

WITH EDGE VORTEX SEPARATION 


Jenn-Louh Pao and C. Edward Lan 


UNIVERSITY OF KANSAS CENTER FOR RESEARCH, INC. 
THE FLIGHT RESEARCH LABORATORY 
Lawrence, Kansas 66045 


NASA Cooperative Agreement NCCI-18 
February 1982 



K'B ikilQB2 

UMGLEY RESEA RPh CENTER 
LIBRARY, NASA 
HAMPTON, VIRGINIA 


l\l/\SA 

National Aeronautics and 
Space Administration 

Langley Research Center 

Hampton, Virginia 23665 



Table of Contents 


page 


1. List of Symbols ill 

2. Introduction 1 

3. Description of Theoretical Formulation 6 

3.1 Mathematical Model of the Wing 6 

3.1.1 Wing Bound Element and Trailing Vortices 6 

3.2 Model of Free Vortex Element System 8 

3.2.1 Leading-Edge Vortex System Geometry 8 

3.2.2 Side-Edge Vortex System Geometry 10 

3.3 Induced Velocity Computations 12 

3.3.1 Induced Velocity Due to Wing 12 

3.3.2 Induced Velocity Due to Free Vortex Element System 14 

3.4 Formulations of Wing Boundary Condition 14 

3.5 Force Free Condition 16 

3.5.1 Formulations of Force Free Condition of Free Vortex 

Element in the Diffused Vortex Filament Model 16 

3.5.2 Formulation of Force Free Conditon of Free Sheet 

Vortex Element in The Core Model 18 

3.5.3 Formulation of Force Free Condition of Vortex Core 

in The Core Model 19 

3.6 Centroid of The Leading-Edge Vortex Filament System 20 

3.7 Z-Coordinates Modification of The Initial Concentrated 

Core 21 

3.8 Formulation of Relaxation Parameters in The Core Model 22 


1 





page 


3.9 Aerodynamic Characteristics 25 

3.9.1 Sectional Aerodynamic Characteristics 26 

3.9.2 Overall Aerodynamic Characteristics 28 

3.10 Summary of Solution Procedure 28 

4. Numerical Results and Discussions 32 

4.1 Numerical Characteristics 32 

4.2 Numerical Results 33 

5. Concluding Remarks 42 

6. References 44 

Appendix A An Improved Lifting Pressure Expression 48 

Appendix B Correction to Sidewash Calculation in Discrete vortex 

Method 56 

Appendix C Determining Initial Core Location by Vector Analysis 66 

Appendix D Correction to Induced Velocities Due to Vortex Core 69 


ii 



w4- w + 


1. List of Symbols 


Symbols 

a 




I 



c 

c 



m 


m 



Percentage of allowed movement of a free segment based on 
the total velocity at Its control point 

Tan a/ Tan (90° - A) 

Fourier coefficients 

Leading-edge boundary condition term 

(x^ - x)t + (y^ - y)j’ + (z^ - z)k m ( ft ) 

(Xj^ - x)i + B(yj^ - y)J + 3(z^ - z)k m ( ft ) 

Induced downwash coefficient due to wing 

Wing span m ( ft ) 

(X 2 - x)i + (y 2 - y)j + (Z 2 - z)^ m ( ft ) 

(X 2 - x)i + 3(y2 - y)l + 3 (z 2 - z)^ m ( ft ) 

Induced downwash coefficient due to leading-edge vortex 
system ( ft“^ ) 

Local chord m ( ft ) 

Mean geometric chord m ( ft ) 

Sectional induced drag coefficient 

Cp - (Cp/Cj^_Q of cambered wings) Induced drag coefficient 

( Drag/(q^S) ) 

Sectional lift coefficient 

Total lift coefficient ( Llft/(q S) ) 

Sectional pitching moment coefficient about Y-axis 
Total pitching moment coefficient about Y-axis based on c 
( Moment /(q_^S c) ) 

The initially computed lift coefficient in the core model 
The reference lift coefficient based on suction analogy 
Pressure coefficient ( )/9_ ) 

CO 00 

iii 



3 + 3 + 


'R 



-i- -J- 

1. J, k 

t 

V 

M 

M 

00 

N 

N 


00 


q 


r 

r 

c 



t 

t’ 



Root chord m ( ft ) 

Sectional leading-edge thrust coefficient 

Total leading-edge thrust coefficient ( Thrust/ (q^S) ) 

Side-edge chord length m ( ft ) 

The Joukowsky force components acting on a segment of 
vortex core N ( lb ) 


a X 


( 


P 


a’ 




»t'l^ IS’ I It'i 

- R) X dt 


)• P 


m ^ ( ft"^ ) 


m ^ ( ft ^ ) 


R. 


Unit vectors along X-, Y- and Z- axes, respectively 
(X 2 - x^)i + (y 2 - yj^)j + (Z 2 - z^)t m ( ft ) 

(X 2 - Xj^)t + S(y 2 - yj^)j + 3 (z 2 - Zj^)k m ( ft ) 
Mach number, or number of spanwise strips plus one 
Free stream Mach number 
Number of chordwlse bound elements 


N(M - 1) 

(M - 1) 

A unit vector normal to the wing surface 

A unit vector normal to the free stream velocity vector 

2 2 2 
Dynamic pressure, (pV^/2) N/m (Ib/ft ) 

2 2 2 
Free stream dynamic pressure (p^ /2) , N/m (Ib/ft ) 

Perpendicular distance of bound element to the core segment 

Characteristic radius of vortex core cm ( in. ) 

p 9 2 2 2 2 ~\h 

[(X2 - + 6 (Y2 - Yj) + 3 (^2 ■ ] 

XX + y^ + z^ m ( ft ) 

XX + 6y^ + 3zic m ( ft ) 

5^ + nj + m ( ft ) 

+ 3nj + 3?^ m ( ft ) 

' iv 



I 


s 



w' 

X, y, z 


» z 

z^ (x, y) 
^min 


2 2 

Reference wing area m ( ft ) 

A unit vector along the free stream velocity vector 

Induced velocity components along X-, Y-, and Z-axes, 
respectively * m/sec ( ft/sec ) 

Free stream velocity m/sec ( ft/sec ) 

Induced velocity vector due to leading-edge vortex system 
} m/sec ( ft/sec ) 

Induced velocity vector due to wing m/sec ( ft/sec ) 

Induced normal wash m/sec ( ft /sec ) 

Wing rectangular coordinate system with positive X-axis 
along axis of symmetry pointing downstream, positive Y-axis 
pointing to right, and positive Z-axis pointing upward, m ( ft ) 

The centroid of the leading-edge vortex filament system m ( ft ) 

The ordinate of camber surface measured from X-Y plane m ( ft ) 

Minimum vertical distance of a vortex segment from the 
wing plane m ( ft ) 


V 



Greek 


a 

0 

Y 

r 



0 

A 

P 

4 ', 

> 1 ^ 


C,n,? 


X 

z 


Angle of attack ( rad ) 

(1 - 

Vortex density referred to freestream velocity 

Concentrated vortex strength based on free stream velocity 
m ( ft ) 

Differential pressure coefficient, 

^lower ^upper 

Chordwlse angular distance ( rad ) 

Sweep angle of wing leading-edge 

3 3 

Fluid density kg/m ( slugs /ft ) 

Spanwlse angular distance ( rad ) 

Local dihedral angle, (Tan ^ ^ 

Sweep angle of leading-edge vortex element or wing 
bound element ( rad ) 

Integration variables in cartesian system m ( ft ) 
Relaxation parameter for the y coordinate of vortex core 
Relaxation parameter for the y coordinate of vortex 


core 



Subscripts 

1 

2 

B 

cp 

i 

j 

k 

Z 

Ze 

L 

R 

t 

T 


The first endpoint of a vortex element 
The second endpoint of a vortex element 
Bound element 
Control point 

Chordwise bound element number 
Spanwise strip number 
Chordwise bound element number 
Leading-edge 

Leading-edge vortex element 
Left trailing leg 
Right trailing leg 
Trailing-edge 

Chordwise trailing vortices 


vii 



2. Introduction 


It is well known that at moderate to high angles of attack, the flow 
separates going around the leading-edge of a sharp-edged wing. The 
separation produces rolled-up vortex sheets above the wing surface which 
emanate from the local leading-edge. The concentrated vorticity above 
the wing induces high suction pressure over a large portion of the upper 
surface of the wing, resulting in a nonlinear contribution to the wing 
loading. In the analysis and design of slender high speed i airplanes 
^^lying at^^f^design conditions, a detailed knowledge of the particular load- 
ing on the wing surface is required to estimate aerodynamic performance 
and structural loads. 

There have been several different theoretical approaches.,in solving 
the vortex flow separation problem since the early 1950 's. In general, 
they can be classified as conical and non-conical flow methods. 

In the first category, based on the slender wing approximation, the 
solutions are obtained from a simplified governing flow equation by 
neglecting the variation of longitudinal velocity. Therefore, the problem 
is solved in the so-called cross-flow plane. Legendre (ref. 1) assumed 
that the rolled-up vortex sheets could be replaced by a pair of concentrated 
cores over the wing and inboard of the leading-edge. Two boundary 
conditions were imposed such that velocities along the leading-edges 
must be finite and the total force on each Isolated vortex must be zero. 

Later, Brown and Michael (ref. 2) used feeding cuts to connect the line 
vortices and the wing leading-edges. Again, the boundary conditions require 
tangency of flow on the wing surface, finite velocity at the leading-edges, 
and zero resultant force on the vortex-plus-cut combination. Mangier and 
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and Smith (ref. 3), modifying the model in reference 2, assumed that the 
vortex sheets are shed tangentially at leading-edges and roll up as spiral 
vortices. The boundary conditions of zero pressure discontinuity across 
the separated vortex sheet and zero total force on the vortex sheet and 
feed cut are imposed. Finally, Smith (ref. 4) developed the model of a 
segmented feeding vortex sheet by improving the numerical procedure 
in reference 3. The main shortcomings of all these models described 
above are limited to conical flow. The assumption of slenderness leads 
to a theory which is Independent of Mach number. Thus, the theory 
breaks down at low Mach numbers in the trailing-edge regions, because 
the tralllng-edge Kutta condition is not satisfied. It breaks down also 
when the aspect ratio is greater than one. 

In the second category, the conical flow assumption is removed so 
that the Kutta conditions are satisfied at both leading-and t railing-edges. 

In the earlier stage of development, the classical lifting surface theory 
was extended by including a relatively simple representation of the leading- 
edge separation. Bollay (ref. 5) suggested a method for low aspect ratio 
rectangular wings in which all the vortlcity is shed from the wing side-edge 
and is assumed to be along straight lines making an angle of a/2 with respect 
to the wing planform. Gersten (ref. 6) and Gamer and Lehrian (ref. 7) 
extended Bollay 's model to cover wings of arbitrary shape. Vorticity is 

assumed to be shed from the leading-edges at an angle of a/2 to the wing 

\ \ 

plane. Th'e wing was replaced by a series of lifting elements. A Multhopp-type 
lifting surface method (ref, 8) was used to solve the Integral equations in 
refierence 7. 

I Nangia and Hancock (ref. 9) combined the conventional lifting surface 
theciry with the detached flow model of Brown and Michael (ref. 2). 
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In the latter model, two concentrated cores lie above the wing, connected to 
the leading-edges by feeding cuts. The trailing wake consists of vortices 
shed from the trailing-edge and the continuing concentrated vortices. To 
determine the separated vortex strengths, their position and the wing vor- 
ticity distributions, the flow tangency condition on the wing surface and the 
trailing-edge Kutta condition are applied at selected collocation 
points, together with the conditions of zero loading on the vortex 
and cut. The leading-edge Kutta condition Is automatically satisfied 
by the choice of the loading function. 

When the wing loading distribution is not required, Polhamus' method of 
leading-edge suction analogy (ref. 10) has been found to be simple and accurate. 
He assumed that the overall lift force can be estimated from the sum of 
potential and vortex lift contributions. The potential lift is defined as 
the lift due to potential flow about the wing with zero leading-edge suction. 
Vortex lift is caused by large suction forces induced by vortices 
in the separated flow, acting on the upper surface of the wing, and 
is assumed to be equal to, in magnitude, the leading-edge suction force in 
the attached potential flow. 

The detailed literature survey of the above flow concepts were 
reviewed by Matoi (ref. 11) and Parker (ref. 12). 

Since 1970, several methods have been developed by including a more 
precise representation of the leading-edge separation. Two popular methods 
under the second category are of the doublet-panel type (refs. 13, lA and 
15) and the free vortices (refs. 16, 17 and 18). In the method of doublet- 
panel type, such as the free-vortex-sheet method of Johnson, et al. (refs. 

13 and 14) , the wing surface and free sheet are replaced with doublet 
panels that have biquadratic strength. The free sheet, whose shape and 
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position must be determined by the iteration, is attached under all conditions 
to a fed sheet. The fed sheet approximately represents the vortex core 
effect and moves with the free sheet. The entire set of doublet strengths 
is determined simultaneously during the Iteration process. Hoeijmakers 
(ref. 15) extended the method of Johnson, et al. (ref. 14) by including the 
entrainment effect of rotational cores. In this model, the spiralling free 
shear layer springing from the leading-edge is replaced by a free vortex sheet 
which is terminated by a combination of a feeding sheet and a discrete line 
vortex and sink. The strength of the sink is related to the entrainment of 
the rotational core and derived from semi-empirical arguments. The combi- 
nation of vortex/sink and feeding sheet sustains no resultant force and pro- 
duces the condition needed for determining the vortex core position in the 
cross-flow plane. From references 13 to 15, the panel method has been shown 
to predict accurate results for flat delta wings. However, 
the computing time involved is lengthy and the methods are generally re- 
stricted to simple geometries. 

In the method of free vortices, Mook and Maddox (ref. 16) have used an 
extension of the vortex lattice method to account for leading-edge separation. 
In addition to standard vortex lattice, a series of kinked vortex lines are 
assumed to be shed from the leading-edge. An iterative procedure is used 
to determine the strength and local flow directions with appropriate boundary 
condition. In the Kandil method (ref. 17), the flow is modelled in a 
manner similar to that of Mook and Maddox (ref. 16) , but included wing-tip 
separation. In this model, the bound elements of the vortex lattices were 
unswept, and the wake behind the tralllng-edge and tip vortex elements 
were force free. The quasi-vortex-lattice potential flow method of Lan 
(ref. 19) has been extended by Mehrotra (ref. 18) to include vortex flow 
effects. This is done by modelling with discrete trailing vortex filaments 
emanating from the leading-edges. The leading-edge Kutta condition is exactly 
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satisfied and partial span vortex separation is allowed. The results of Meh- 
rotra's model (ref. 18) indicate that the predicted spanwise pressure distri- 
bution is more diffused than the experimental data show. This is because the 
effect of a concentrated vortex core in the real flow is not well represented 
by a number of free vortices. From the published results, it is clear that 
none of the existing free-vortex methods is capable of predicting accurately 
the lifting pressure distribution over the wing surface within a reasonable 
amount of computer time. 

The main objectives of this investigation are : (1) to reduce the computer 

time through the use of a pair of concentrated vortex cores above the wing sur- 
face, which should also improve the pressure prediction of Mehrotra's model 
(ref. 18); (2) to allow the effect of side-edge vortex separation; and (3) to 
extend the method of reference 18 to more realistic geometries. Configurations, 
which are applicable to real airplane geometries include features such as 
camber, strakes and leading-edge vortex flaps. All these configurations will 
be treated in this report. 

In the following, the aforementioned improvement will be described. First, 
the leading-edge vortex system is superimposed on the regular quasi-vortex- 
lattlce grid by using segmented vortex filaments, as shown in reference 18. 
Second, the side-edge vortex system is obtained by extending all the bound 
elements to the tip chord in the last vortex strip. Third, the leading-edge 
diffused vortex elements, in reference 18, are allowed to merge into a con- 
centrated core after two iterations. 
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3. Description of Theoretical Formulation 


Mehrotra's method (ref. 18) is a diffused vortex-filament model of the 
vortex core but serves as a basis for the present method. In the present 
analysis, the basic assumptions are: (1) the wing is represented by a bound 

vortex sheet across which there exists a pressure difference; (2) the sep- 
arated flow along leading-edges is represented by force free leading-edge 
vortex elements which feed vorticity into a concentrated core. In the present 
method, the Quasi-Vortex-Lattice method (QVLM) (ref. 19) is used to simplify 
the induced velocity expressions due to the wing bound vortex sheet. 

During the iteration process, the following boundary conditions are 
imposed: (1) the wing surface must be Impermeable; (2) Kutta conditions are 

imposed along the leading- and trailing-edges of the wing; (3) in the diffused 
vortex-filament model, the force free condition is applied on the leading-edge 
vortex filaments and trailing wake elements; and (4) in the core model, the 
force free condition is applied on the free vortex elements , concentrated 
core and trailing wake elements. 

This is a non-linear problem because the strength of the wing bound vor- 
tices and free element vortices, the locations of the leading-edge free 
vortex elements, and the concentrated core are unknown. 

3.1 Mathematical Model of the Wing 

The wing lies in the x-y plane and the x-axis is taken to be streamwise 
along the wing center line. The origin of the right-handed rectangular coor- 
dinate system is assumed to be at the moment reference point of the wing. 

3.1.1 Wing Bound Element and Trailing Vortices 

According to the QVLM, the semi-circle method is used to determine the 
control station locations and the vortex strips in the spanwise direction, and 
the wing control points and bound elements in the chordwise direction. The 
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locations of bound- and tralllng-vortex elements are as shown in Fig. 1. The 
x-location of bound elements is given by 


x^ = Xj^ + -| (1 - cos (— 2 n '~ '^)) ♦ ^ N 

where x^ is the leading-edge x-coordinate, c is the local chord and N is 
the number of bound elements in the chordwise direction. The spanwise lo- 
cations of trailing vortices are given by 


rl 3 - - .J - ^))^ j = 1, M (2) 

where b is the wing span and M is the number of legs of trailing vorticity, 


Yj ■ 4 <1 - <- 2M 


which is one higher than the number of spanwise vortex strips. The locations 


of control points are given by 

k 

x__ = x„ + (1 - cos k = 0, 1, N 




-4^ '2 


^cp " j =■ 1. (M - 1) 

where x- and c. are the leading-edge x-coordinate and local chord at y , 
3 ^ ‘'Pj 


(3) 

(4) 


respectively. 

For the strake-wlng configuration, there are two sections in the span- 
wise direction. Each spanwise section is divided into vortex strips by the 
semi-circle method. The edges of the vortex strips in each section are ob- 
tained through the following relation 


yj = y'k-i + 2^ (1 - j — \ 

and y control points are given by 

Ycp^ ' Yfc-1 +r ^ 

where b^ is the width of the k^^ spanwise section and is the number of 
spanwise strips plus one in k^^ section (see Fig. 2). 


(5) 


( 6 ) 
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t 

y, 1 are defined as 
•^k-1 


I 



3.2 Model of Free Vortex Element System 

3.ZJ. Leading-Edge Vortex System Geometry 

The leading-edge vortex system is superimposed on the regular quasi- 
vortex-lattice grid. A typical leading-edge vortex element is shown by 
points A through J in Fig. 3. These points are connected by a series of 
short straight segments. A typical set of Initial locations for these 
segments is indicated with dashed lines and the locations after two iterations 
are given by solid lines. 

The segments of the vortex element have the following characteristics: 

(1) Initially, points A through E lie along a wing trailing vortex 
filament with point A being one root chord away from the 
trailing-edge (point D) in the downstream direction, and the line 
segments between A and D are parallel to the axis of symmetry. The 
segments B-C and C-D are 0.1 c^^ long. The length of each line segment 
between A and B remains unchanged as the solution progresses. Segment 
B-C is allowed to move only in the vertical direction, whereas segment 
C-D is fixed in the wing plane because of the trailing-edge Kutta 
condition. Segment D- E is also fixed in the wing plane. 

(2) Points E, F, G and H lie in the wing plane. The location of segment E-F 
is ahead of the wing first bound element and is given by 
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X 

E 


(7) 


X 


y ( 1 - COS ( 


2 (N+1) ^ ^ 


F £p 2 ^ 2 (N+1) 


) ) 


( 8 ) 


The segments F-G and G-H are of the same length and point G lies on the leading- 
edge. These segments are fixed in the wing plane to satisfy the leading-edge Kutta 
condition. 

(3) The initial location of point I is a function of a and is given by 
(ref. 18) 

(9) 

( 10 ) 


X = X 

I F 




Zj. = 0.1 Cj^ Tan (22.5 - 0.5 a), for 15 


(J.1) 


or 


= 0.1 Tan a, for a > 15° (12) 

V/here Cj^ is the root chord and a is the angle of attack. Initially point 
J is one root chord away from the tralllng-edge. The line connecting 
I and J may be divided into segments of equal length which lie in a 
plane parallel to the x-z plane. The initial height of these segments is 
approximately 0.1 c^^ above the wing plane. In the final converged location, 
all segments between point H and J (free vortex elements, concentrated core 
and trailing wake elements) are aligned in the local velocity direction at 
their mid-points to satisfy the force free condition. 

(4) The semi- infinite segments from point A to infinity and J to infinity are 
straight and are parallel to the undisturbed free stream direction. Since 
the far wake of the leading-edge vortex element has very little effect 
on the wing aerodynamic characteristics, the far wake is assumed to be semi- 
infinite and parallel to the free stream. 



3.2.2 Side Edge Vortex System Geometry 

To establish the side-edge vortex system as shown in Fig. 4, the leading- 
edge vortex element and wing bound element of the last vortex strip is modified 
as follows: 

(1) Extend the leading-edge bound element to the side-edge (point D) . Points 
D, E, F, G and H lie in the wing plane. The segments F-G, G-H, H-I and 
I through J are given in the same manner as in the preceding section. The 
extended leading-edge bound element consists of segments D-E and E-F. 

The location of each point is given by 

■ %ip ( 1 - cos ( ) (13) 


-F - + -r <1 - < TW)> '> 

r 


^E 2 ^F^ 


(lA) 

(15) 


(2) The initial location of each vortex segment along the trailing leg of 


the extended leading-edge bound element, point A through D, is given by 


Xj^^ = x^ + Ay cos 2 


(16) 


^j+1 ^Tip 


"j+r + Ay sin 2 


(17) 

(18) 


where 


Ay = 0.95 X (width of the last vortex strip). 

(3) Extend the regular wing bound element to the side-edge. (Points P, Q, R, 
S) . Points P, Q, R and S lie in the wing plane and at the side-edge. The 
location of each point is given by 

X. = X 


X. 

Txp 


^i = Xlp 


. *^ip /■, ,2i - 1 

+ — (1 - cos (■ 


2N 


17 )), 1 = 1, N 


(19) 

( 20 ) 
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where x. is the leading-edge x coordinate at the side-edge of the wing; 

\ip 

the side-edge local chord length. 

(4) Initially, all the segments between points D and A, P and T, Q and V, etc. 
lie in a plane parallel to x-z plane at the wing side-edge. Points A, T, 
U, V and W are one root chord away from the trailing-edge. The initial 
side-edge vortex segment along points P through T (or Q through U, etc.) 
is given by, 

= x^ + Ax cos y (21) 


^1+1 ^Tip 


( 22 ) 


z^^j^ = z^ + Ax sin Y (23) 

where Ax is assumed to be half length of the side-edge chord. 

(5) The semi- infinite segments from points A, J, T, U, etc. to infinity are 
straight and are parallel to the undisturbed free stream. 
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3.3 Induced Velocity Computation 
3.3.1 Induced Velocity Due to Wing 

In the quasi-vortex-lattice method (ref. 19), to satisfy wing boundary con- 
dition, the continuous vortex distribution over the wing is replaced by a quasi- 
continuous one, being continuous chordwise but stepwise constant in the spanwise 
direction. Thus, the wing surface can be divided into a number of vortex 
strips with the associated trailing vortices (Fig. 1). In any strip, consider 

a vortex element y dx with an arbitrary direction i (Fig. 5). The induced 

til 

velocity due to all bound elements in the i strip is given by (ref. 19) 


ll 4TT X^ 


') 


a X t 


( 




_ .t'dx’ (24) 


fg'l P'l 

and due to the associated trailing vortices by 

(R) = 3^ f^t Y (^') (^i - ^) X dt ) 

^2 X. x’ „3 

^ (25) 

The transformation, x* = + c(y) (l-cos0 )/2, reduce equations (24) and 

(25) to 

V (R) = 

^1 


and 


where 


(4>) = 


8it 

0 

6^c(y) 

./ 

8 •n 

0 

a X 



B^c(y) f’' ( 6) Y ( e) SinG d 0 

'^i, ■ “O 0 


^ f 

( )• 


|a-xt'|2 |^,| |-,| 


s, (♦) -JT 


(26) 


(27) 


and c(y) = 3 ^ - x a, etc. are defined in the List of Symbols. 
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,th 


The total induced velocity due to the i strip of vortex distribution 
is given by 

2 TT 

Vi (h = J (e) Y (e) Sine 


de 




de 


- g c . (y) / (e) y(e) sine de 

OTT 0 Z 


(28) 


Where the first term is due to bound element, second due to left leg of trailing 
element and third due to the right leg of trailing elements. The above integrals 
are reduced to finite sums through the mid-point trapezoidal rule (ref. 19): 


N 


^1 *2^ - % > Si" 


where = (2k - 1) it/ 2N and locations of bound elements are given by 






2V - 1 

= h a - Cos ( ~ 2^ — TT)), k = 1, N 


and 


x» _ the leading-edge X -coordinate at y, (left leg) 
1" ^ 

"x. ^ = the leading-edge »-cordinate at (right leg) 


(29) 


(30) 

(31) 

(32) 


Cj^ = chord length at y^ 

= chord length at y^ . 

The location of control points and trailing vortices are described in equations 
(1) through (6). 

Thus, the Induced velocity due to all vortex strips of the wing can be 


written as 



3.3.2 Induced Velocity Due to Free Vortex Element System 

The leading-edge free vortex system consists of M - 1 distinct vortex 

elements. Each element may have different number of small vortex segments, 
til 

Assume that the i set has L small segments. The induced velocity at a given 

til til 

point (x,y,z) due to j segment of i element is given by (ref. 19) 




t 




a'xt’l^ \P 


-) .t' 


(34) 


where 


^ 

R = xi + y j+zk 


a = (Xj-x) T + (y^ - y) J + (z. - z) k 
^ = (Xj+i“ x) 1 + (Yj+i- y) J + (2j+i“ ^ 


^ = (Xj+i - ^j> ^ + (yj+r y) 3 + z.) t 

a' = (x^ - x) 1 + 3 (Yj - y) 1 + 3 (Zj - z) S 
1^’ = - x) 1 + 3 (Yj+i" y) 1 P ^^j+1 " ^ 

P = (x^+i - X.) t + 3 (y.^i - y.)!” + 3 z.) t 

til 

= leading-edge bound vortex strength of i strip. 

The subscripts j and j + 1 correspond to the end point of j‘‘“ segment. The 

induced velocity due to i^^ element can be written as 

L 

0. (J) =y~ S , , (5) (35) 

1 fir 

Therefore, the induced velocity due to the entire leading-edge free vortex 
elements is 

U (t) = U,. (^) (36) 


.th 


i=l 


3.4 Formulations of Wing Boundary Condition 

The bound elements and the corresponding control points of the wing 
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surface are numbered from the leading-edge to the traillng-edge and from 


the root to the side-edge. The flow tangency condition can be written as 


‘V'i 


dz 


Cosa ^ ^ i “ Sina 




N 


4 - 


N, 


( 37 ) 


where is the total number of wing bound elements; the number of 

leading-edge vortex elements; the induced downwash at the i^^ control point 

of wing due to a unit horsehoe vortex density at j ; the induced downwash at 

the i control point due to the k leading-edge vortex element of unit 

th 

strength; y. the vortex density of the j bound element; P, the strength of 
2 

til th 

the k leading-edge vortex element; (dz/dx)^ the camber slope at the i 

control point and a the angle of attack. According to equations (52) and (53) 

of reference 19, the leading-edge thrust coefficient is related to the normal- 

wash on the leading-edge by the following equation 


dz I 

a , = (induced normalwash - Cos a (-;^). + Sin ot) | 

JvcC QX k 3t • S • 


(38) 


which leads to the last expression on the right hand side of equation (37). In 

the expression, is defined as (ref. 19) 

2 c Cos A 
k 


£k 


= N J Tan^ A + { 


/ 2 

TT Jl - M Cos A 
“ 00 


where 

N = number of chordwise bounded vortex elements 

A = leading-edge sweep angle 

M = free stream Mach number 
00 

e =/i - 

Y 00 


th 


= sectional leading-edge thrust coefficient of the k strip 


For complete leading-edge separation cases, the sectional leading-edge suction 
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Is zero and so is a 


Jlk* 

3.5 Force Free Condition 

The vortex segments of each free element above the wing surface and the 

trailing wake are to be aligned in the direction of local velocity vector cal- 

til 

culated at their mid-points. Consider the i segment of a vortex element. 

The coordinates of its end points are given by (x , , y , , z.) and (x. , , , Yj » 

. Assume that the induced velocity at the mid-point of this segment at 
a given iteration is given by (ui + vj + \vk). In the following sections, the 
formulations of force free condition will be derived in the diffused vortex 
filament and core models respectively. 


3.5.1 Formulation of Force Free Condition of Free Vortex Elements in The 

Diffused Vortex Filament Model 

til til 

For the i segment, the new location of the (1+1) end point can be ob- 
tained by aligning the segment in the direction of local velocity vector. 

That is 


^ 1+1 

^i+l 

"l+l 


= X. + -^ As 
X U 

1 V . 

= Yi + U As 

. w . 

= z^ + y As 


(40) 

(41) 

(42) 


where 

U = (u^ + v^ + w^)^ 

*= ■ + <^ 1+1 - 
Before equations (40) - (42) are used, the following points should be considered: 

(1) The length of each segment is to be preserved. 

(2) The free vortex segments above the wing should not come too close to 
the wing surface to avoid numerical difficulty in the present inviscld 
theory. 

(3) The adjustment of the location of each segment to satisfy the force 
free condition should be such that it does not result in numerical 
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fluctuations. 


Based on the above considerations, equations (40) to (42) will be modified as 
follows : 

Consider the same i^^ segment. If this segment moves "a" percent only 
according to the velocity computed at its mid-point, the equations (40) to (42) 
can be modified to be 


Ay' = 

V 

^ u 

As 

+ 

(1 - 

a)(yi+i - 

-yf) 

(43) 

Az' = 

w 
^ U 

As 

+ 

(1 - 

a) (z^^^ - 

- z.) 

(44) 

Ax' = 

|as' 

1 

(Ay')^ 

- (&z')2 


(45) 

i that 








^1+1 " 

' X, 
X 

+ 

Ax' 




(46) 

^1+1 = 

■ ^i 

+ 

Ay' 




(47) 

^i+1 = 

' "i 

+ 

Az' 




(48) 


The free elements of the leading-edge vortex system have been restricted not to 
come any closer than a minimum specified height to the wing surface, which is 
given empirically by Mehrotra (ref. 18) as follows 

(49) 


z . = 0.1 c_ Tan (22.5 - 0.5a) , for a < 15 

min R — 


z = 0.1 c_ Tan a 
min K 


, for a > 15 


(50) 


If z . . , is less than z . , it is then set z,,, equal to z . and Az' is re- 
i+1 min 1+1 ^ min 

calculated by using 

■ ^mln - "l 

This restriction is needed because whenever the free elements are close to the 
wing surface, they induce large velocities on the wing and vice versa, which 
makes the locations of free elements fluctuate. In the real flow, at small 
angles of attack, the leading-edge vortex system is weak and diffused. The 
effect of the free vortices is artificially reduced by increasing as the 

angle of attack is decreased below 15 degrees. Note that equations (43) through 
(48) described above are mainly used to (1) establish the leading-edge vortex 
filament system in the first two iterations; and (2) find the new position of 
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trailing wake elements and side-edge vortex system in the entire iteration 
process. 


3.5.2 Formulation of Force Free Condition of Free Sheet Vortex Element in The 


Core Model 
til 

For the i segment along a free vortex element between the leading-edge 
and concentrated core, the initial slopes are determined by the total in- 
duced velocities at the mid-point of this segment 

(52) 


'•Ax'^ u 


.Az . ' _ w 
''Ax^ " u 


• th 


The new slopes of the i segment of a free vortex element are given by; 

(^)* . d* (^) + (1 - d*)(^)' 

Ax Ax Ax 


K7 ^ 'k Ky ^ 


(53) 

(54) 

(55) 


Similar to the parameter "a" in equations (43) and (44) of section 3.5.1, the 

k 

relaxation parameter of free sheet, d^ is used in equations (54) and (55) in 

the core model, to avoid any rapid variation of free sheet shape and divergence of 

* 

the iteration procedure. The initial value of d is 0.5. During the itera- 
tion process, when the rate of change of the total free sheet vortex element 

* 


force is less than five percent, d is set to 0.95. 


Since the length of each segment is conserved, therefore, 

,Az.2 -\h 


.3 - A. [1 ^ 


(56) 


Hence 


* 

(Ax) 


It follows that 


As 


(Ay)* = (Ax)* (||)* 
(Az)* = (Ax)* (f)* 


(57) 


(58) 

(59) 
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Thus, the new location of the (i+l)*"^ end point is given by 

* 


+ (Ax) 

^i+l = ^1 + 

"i+1 = "i + 


( 60 ) 

(61) 

(62) 


3.5.3 Formulation of Force Free Condition of Vortex Core in The Core Model 
During each Iteration, the Joukowsky force acting on each segment of the 
core is computed. Let 


F = V X r 

= (ui + ^^ + w^) X ^ 

Then 

F r A 

y _ w z _ _w ,Az. 

r u u r u “ ^Ax-' 

X X 

F r . 

— = _V 

r u r u '■Ax'' u 

X X 

If the new position of the core is such that F =0 and F =0, then 

O' - O' ■ . ’ ■ 

(My’ _ (V)' = 0 
'■Ax'' 'u'' 


(64) 

(65) 

(66) 
(67) 


Assume that 


0 = 0 
= O’ 


Combining equations (64), (65) and (66), (67) yields the following equations 


for the l*"^ segment : 


(M)’ = (M) - (__Z) 

'Ax^i ''Ax-'i 


( 68 ) 


(— ) ' = (— ) +( — ^) 
'■Ax''i ^Ax''! u% 

± -L 3 J 1 

,th 


The new slopes of the i segment of the concentrated core are taken to be: 

(th = (1 - \ €>i 

Where Xy and are relaxation parameters and are discussed in Section 3.8. 


(69) 

(70) 

(71) 
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Combining equations (68) through (71) results in: 


X X XX 

O! - 

X 1 X 1 


( 70 ) 


(71) 


Daring the iteration process, when the rate of change of the vortex core force 
is less than five percent, the relaxation parameters X and X„ of each core seg- 

y ^ 

ment will be further reduced by multiplying the factor 0.1. Thus, the move- 
ments of core segment will not cause excessive fluctuation in the aerodynamic 
characteristics . 


Since the segment length is conserved, it follows that 


As = Ax 


+ <S>' 


]' 


(Ax)* = 


As 




(74) 

(75) 


Hence, 


(Ay)* = (Ax)* (|^)’ 


(Az)* = (Ax)* (^)* 
i 


(76) 

(77) 


,th 


Thus, the new position of (1+1) end point of the concentrated core are found 
as follows 

■k 

(78) 

(79) 

(80) 

In the present analysis, the effect of the secondary vortices has not been in- 
cluded. 


^1+1 ^ ^i 
5'i+l ■ >'i + 
"l+l " "l 


3.6 Centroid of the Leading-Edge Vortex Filament System 

After establishing the leading-edge vortex filament system based on exer- 
cising Mehrotra's code (ref. 18) through two iterations, the centroid of that 
diffused vortex element system can be found as follows (ref. 20): 
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Consider a series of cross flow planes, proceeding from the wing apex 

toward the trailing-edge. The centroids of the vortex filaments penetrating 

these planes are given by, 
n . 



z . 

3 


iL ^i ^ij 

r . 

nj 

s r z 

i=l ^i ij 


cj 


(81) 

(82) 


where , 


, Zj are the centroid position in the j plane, 

til 

Yy , z^j are the intersection position of the i vortex filament with the 


jth plane. 


th 


is the circulation around the i vortex filament, 

th 

r . is the concentrated line vortex strength at the j plane, 
n^ is the number of vortex filaments penetrating the j plane. 


3.7 Modification of Z-Coordinates of the Initial Concentrated Core 

In order to have better starting solution so that the number of iterations 
and computer time can be reduced, the initial z-coordinates of the concentrated 
core from Section 3.6 are modified to take the experimental values summarized 
by Smith (ref. 4): 

z = 0.154a' + 0.1 a' < 1.2 

V 

z^ = 0.1333a' + 0.158 1.2 <_ a' < 2.2 (83) 

z = 0.13a' + 0.132 a' > 2.2 

V 

where 

a' = Tan a/Tan G 
0 = 90° - A 

z^ is the nondimens ional z-coordinate of the concentrated core, referring to 
the local semi-span; A the leading-edge sweep angle; and a the angle of attack. 
Equation (83) will be used only in the simple wing configuration. 
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The initial configuration oh the core model, on the right half of the 
wing, is shown in Figs. 6A through 6C. The top view of the Initial free sheet 
and the modified y-coordlnate of the core location are shown on the x-y plane 
in Fig. 6A. The initial modified y-coordinate of the core location is such that 
it is at least at a location equal to 60 percent of the local semi-span of the 
wing (see references 4 and 22). The rear view of the detailed free sheet shape 
and its relative vortex core location above the wing surface are shown on the 
y-z plane in Fig. 6B. The side view of the free sheet and the modified z-coor- 
dinate of the core location are shown on the x-z plane in Fig. 6C. Note that 
the initial modified z-coordinate of the core location is based on equation (83). 


3.8 Formulation of Relaxation Parameters in the Core Model 

In the first two iterations of the core model, the core is forced to 

move either in the y or z direction. The direction depends on the sign and 

the difference between the initially computed lift coefficient and the reference 

lift coefficient based on the suction analogy (references 26 and 27) . From 

the numerical experimentation, the relaxation parameters X and X , in equa- 

y 2 

tions (70) through (73), are taken as follows; 

First iteration 

If the initially computed lift coefficient is smaller than the reference 
lift coefficient based on the suction analogy, the core will move down and 


the corresponding rel^gatlon^pajameter 

li Xj • 

X =-15 s ^ 
z, (C. Tan a) 

^s 


If the initially computed lift coefficient is greater than the reference 


lift coefficient, the core will move out and the relaxation parameter is 


(Cj^ ~ ^ ° 


X = -2.0 


where is the reference lift coefficient, and is the initially computed 
^s 1 

lift coefficient. The initial set of X and X are designed so that the 

^1 ^1 
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movement of core segments is proportional to the difference in (L and . 

^i s 

Also, the higher the angle of attack, the larger the value of X or the lower 

^1 

the value of X is used. It was determined from numerical experimentation that 
^1 

the initial core needs larger move-down movement in the z direction than the 
move-out movement in the y direction. 

Second iteration 

If the initial movement of the concentrated core is in the z direction, the 

y-coordlnates of the core will be adjusted next in the second Iteration. The 

computed lift coefficient after the initial adjustment of the concentrated core 

is denoted by (^. The core will be caused to move out if is greater than 

and to move in if is less than Cj^ . The corresponding relaxation para- 
s ^s 

meter is: 


X = - 


^s 


( 86 ) 


On the other hand, if the initial movement of the concentrated core is in 
the y direction, the z-coordinate of the core will be adjusted in the second 
iteration. The core will move down if is less than and move up if 

L Jj_ Li 

S 

is greater than . The corresponding relaxation parameter is: 


(C^ Tan a) 


(87) 


The relaxation parameters X and X have the same form as the initial set ex- 

cept the magnitude is reduced by half. Special conditions are imposed for all 

configurations when the angle of attack is below 20 degrees. They are (1) 

the magnitude of X in equation (84) is to be reduced by half; and (2) for the 
^1 

angle of attack below 10 degrees, the angle of attack in equations (84) and (87) 
is set to 10 degrees to avoid the singularity in a in computing X and X 

Zi 

Therefore, the movement of the core in the first two iterations for angles of 
attack below 20 degrees will not cause a large fluctuation in the aerodynamic 
characteristics . 
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The relaxation parameters X and X of each segment of the core are to be 

y z 

used from the third iteration on, based on the force gradients of core segments 
before and after the adjustment of the free sheet vortex elements. Thus, the 
relaxation parameters to the i*"^ segment of the core are: 


X = C ~ 

>'1 ^1 'V' i 

F - F 

- <=Z (If > 

z . z . r . 

1 1 zl 1 




1= EF 


( 88 ) 


(89) 


In the above, ^^l’ ^y2 ^z2 segment forces of the core before 

and after the adjustment of the free sheet vortex elements, and EF is the 

algebraic sum of segment forces acting on the core after adjustment of free sheet 

vortex elements. The sign of C and C are determined so as to reduce the 

“ y . z . 

•'i 1 

total force acting on the free sheet vortex elements. 

In the core model, the locations of free sheet vortex elements and concen- 
trated core have a very sensitive effect on wing loading. To avoid the div- 
ergence of the iteration procedure: (1) the restrictions are imposed on the 

movement of free sheet vortex elements and concentrated core as discussed in 
sections 3.5.2 and 3.5.3; (2) the end point of the first core segment near 
apex will remain fixed from the third iteration; and (3) the z-coordlnates of 
concentrated core are not allowed to be below given by equations (40 and (50). 


24 



3.9 Aerodynamic Characteristics 

In reference 26 , the more accurate expressions to calculate the overall 
aerodynamic characteristics are derived by Lan for the cambered wing. Both 


9z 


and 


3z 
c 

3y 


3x 


are included in the computation. 


Assume that the wing surface (Fig. 7) can be described by: 
z = z^(x,y) 

where z^(x,y) is the ordinate of camber surface measured from the x-y plane. 
Introduce a function f such that: 
f = z - z^(x,y) 

Therefore, a unit normal vector on the wing surface can be defined as: 


(90) 


(91) 


n 


Vf 

|^f| 


3z 3z 

c -f c -f . t 

3x ^ ~ 3y ^ ^ ^ 


3z 2 


3z 2 nk 


+<3y 


dz z n 

] 


(92) 


The unit vectors n^ and t^ are related with n (Fig. 8) in the following expressions 
n^ = -Sin a 1 + Cos a k (93) 

= Cos a ? + Sin a ^ (94) 

The computed lifting pressure AC^ (see Appendix A) is acting normal to the local 
camber surface, i.e., in the n direction. Then, this pressure force can be de- 
composed into n and t directions to determine the lift and drag components. 


(95) 


n*t = ( - 


h 


respectively. Hence, the component in the lift direction is: 

3z r 3z 2 2-j H 

the component in the drag direction is: 

3z — 3z n 3z ft *1 

3 ^ Cos cx + Sin a)/ [l + + (^) J 

To find the pitching moment about y axis, the j component of (Ax i + Az k) x n 
is needed, where Ax and Az are the moment arms of the pressure elements refer- 
enced to the coordinate system. In the present analysis, Az is assumed to be 


zero. Hence, the corresponding j-component of (Ax i x n) is: 

c 2-lJs 


r 2 2n 


(97) 
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3.9.1 Sectional Aerodynamic Characteristics 

To find the sectional characteristics, it is assumed that the wing vor- 
tices are situated along the cambered surface. Since the resulting lifting 
pressure is acting normal to the local camber surface, the sectional lift. 
Induced drag and pitching moment coefficients can be determined by combining 
the lifting pressure AC^ and equations (95) to (97) through the local chord 
Integration. Hence, the sectional lift, induced drag and pitching moment co- 


th 

efficients at the j strip have the following expression: 



r 

3z _ 3z _ .. 


AC (t — Sin a + Cos a)/ 

3 P 3x ' 1 

.1+ <3y ) ] 


dx 


[ ■ 9z ^ P 3z fy 3z nq 


(98) 


3z 

AC ( - -^ 
p 3x 


3 

P 3z 

, [- 

-i “ “I 


Cos a + Sin 


- c 


Cos a + Sin 


a)/ [ 


3z 


3z 


1 


dx 


3z 




2 ^^c 2 

j ^ 


(99) 


c 

m 


‘j 


cX r 

3z ry 9*1- 

-J AC - (Ax)/ 1 

■^x„ ^ 

1 + <3. ) + ] 


(100) 


where the terns in equations (98) and (99) are the leading-edge suction 
3 

effect; Xjj^ and x^. the leading-and trailing-edges x-coordinates of the chord 

J th _ 

passing through the control station of the j strip; c the mean geometric chord, 

Cj the local chord length. The subscript in equations (98) and (99), means 
3z 3z^ 

(-g^) and (gy^ ) are evaluated at the leading-edge of the j strip. Using 
3 3 

this transformation. 


X 


= Xj^ + (1 - Cos 6) 

3 


equation (98) can be reduced to: 

3z 


= "'Sj ’'ACp Sin 0 (-^ Sin a + Cos a)/[l + (-^) + (-^) ] 


d0 
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+ C 


(- 9z 


) Sin a + Cos a J j ^ 


9z 


I + + (^)^1 


j" e 
9z o 9z 


N+1 p 9z P 9z « 9z j-,h 

7^^Tt 2I AC (-T-^) Sin a + Cos a Sin 0, / 1 + (-^) + (^) 

(N+1) Pk L 3x J k' L '9x '9y ^ J 


[ 3z ^ p 3z « 3z *-• 

(^) Sin a -f Cos a]^ / [i + ^ + (^)2] 


h 


3 e 


(101) 


and 


0 (2k - 1)TT 

k 2(N + 1) 


k=l, (N + 1) 


where the mid-point trapezoidal rule has been used to reduce the integral into 

a finite sum (ref. 19). Similarly, the sectional induced drag and pitching 

.th 


moment coefficients for the j strip are given by 


N+1 p oz P oz 

dj - I- “J V f + <a;r> 

3z r T r 9 

■ "t. [- 

J'k-' 3^ J -lei J 

] ' 

-‘z 

e 

N+1 c. . az 

V]-Vh^: 

? 


9z , 

''9y ■' 


( 102 ) 




(!!c)21% 

^3y 


(103) 


According to equations (52) and (53) of reference 19, the sectional 
leading-edge thrust coefficient is given by: 


2 2 fw! -(Cos a (-^) , - Sin a)"! 

_ u(l - Cos A) L j ^ dx 3 ^ 


“j 


2N^ Cos A (1 - m£ + Tan^ A) 


(104) 


where is the free stream Mach number; A the sweep angle of the leading-edge; 
dz 

wl and (-T— ) are the induced normalwash and slope of the wing surface at the 
3 dx J 

leading-edge. For sharp-edged wings, is assumed zero. 
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3.9.2 Overall Aerodynamic Characteristics 

The total Induced drag, lift, pitching moment and leading-edge thrust 
coefficients of the wing are determined by spanwise integration of the sectional 
characteristics. Again, the integration is first transformed to an angular 
coordinate (0 ^ ^ tt) , and then reduced to finite sum by the conventional 

trapezoidal rule (ref. 19). Therefore, the total lift coefficient has the 
following expression, 

'l ■ I V (105) 

Jo 


where b is the wing span and S the wing area. By the transformation: 
y = (1 - Cos <j>) 

equation (105) can be reduced to: 

b 

S. = 2S J Sin <1. d<{. 

o 


bir 

~ 2SM 


M-1 

c c. Sin (b. 


and 

<l>j^ = , i=l, (M-1) 

where M-1 is the total number of spanwise strips. 

Similarly, the Induced drag, pitching moment and leading-edge thrust 
coefficients are given by: 


(106) 


(107) 


s." 

M-1 

Hi c, sin 

(108) 

1 

X=1 1 


c s 

m 

2^ ^ '^m . ^^i *^1 

(109) 

r — 

Ct - 

i=l i 
M-1 

^ . ^i *^1 

1-1 1 

(110) 


3.10 Summary of Solution Procedures 

In the diffused vortex filament model of Mehrotra (ref. 18), the basic 
unknowns of the problem are the bound vortex density on the wing, the strengths 
and locations of the elements of the leading-edge vortex system and the 
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trailing wake. The problem is solved in the following iterative manner: 

(a) Prescribe the vortex-lattice for the wing surface, and the 
initial locations of the leading-edge and side-edge free 
elements over the wing and in the wake. 

(b) By satisfying the wing boundary condition, equation (37), ob- 
tain the bound vortex density of the wing and the strengths of 
free elements. 

(c) Calculate the aerodynamic characteristics and loads. 

(d) Adjust the free elements of the leading- and side-edges vortex 
system and the trailing wake in the local velocity vector dir- 
ection. 

(e) Repeat steps (b) through (d) for two iterations. 

After two iterations in Mehrotra's diffused vortex filament model, the 
present model with the leading-edge vortex core is then introduced (See Fig. 9), 
In the core model, the basic unknowns of the problem are the bound vortex 
density on the wing, the strengths and locations of the free sheet vortex 
elements, the concentrated core and the trailing wake. The problem is still 
solved in the iterative manner: 

(f) Find the centroid of the established leading-edge vortex 
filament system, equations (81) and (82), from step (e) . 

(g) Allow the leading-edge vortex elements to merge and feed 
vorticity into the concentrated core through connecting 
segments. 

(h) From the numerical experimentation, the z-coordinate of the 
initial core is modified by equation (83), and the y-coordlnate 
of the initial core is adjusted so that the initial core 
location is at least 60 percent of the local semi-span 

of the wing to have a better starting solution. 
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(i) By satisfying the wing boundary condition, equation (37), 
the initial bound vortex density of the wing and the 
strength of free vortex elements are obtained. 

(j) In the first two iterations of core model, the core is 
forced to move either in the y or z direction. The dir- 
ection depends upon the sign and the difference between 
the initially computed lift coefficient and the reference 
lift coefficient based on the suction analogy. The relax- 
ation parameters and X^ are based on equations (84) 
through (87). 

(k) From the third iteration of the core model, the movement of 
the core is based on the relaxation parameters for each seg- 
ment of the core, equations (88) and (89). These are used 
to keep the core moving in a direction for which the total 
algebraic sum of the forces acting on the free sheet is to 
be decreased. 

(l) Based on step (k) , the orientation of the leading-edge free 
elements of the free vortex sheet, side-edge vortices, the 
concentrated core and the trailing wake are adjusted in the 
local velocity vector direction. 

(m) By satisfying the wing boundary condition, equation (37), 
the bound vortex density of the wing and the strength of 
free elements are obtained. 

(n) Calculate all the aerodynamic characteristics and loads. 

(o) Calculate the force acting on all free elements of the free 
Vortex sheet over the wing surface and the concentrated core. 

(p) Repeat steps (j) through (o) until a converged solution is 
obtained. 
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The converged positions and orientations for the free sheet vortex 
elements, side-edge vortex elements, and the concentrated core are assumed 
to have occurred when the sum of the magnitude of the total force acting on 
the core and free sheet has reached a local minimum value in the iteration 
cycle. 

Numerical results to date indicate that when the angle of attack is 
decreased below 20 degrees, the present method has difficulty in generating 
a satisfactory starting solution. The main reason is the poor initial shape 
of the free vortex sheet from Mehrotra’s model. Thus, with the core model 
method, the starting solution for wings at angles of attack less than 20 de- 
grees is taken to be the initial shape for 20 degrees angle of attack. 
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4. Numerical Results and Discussions 


This section presents numerical characteristics of and predictions by the 
present method. These predictions are compared with results from the methods 
of Mehrotra (ref. 18), Johnson, et al. (refs. 13 and 14) (called Free Vortex 
Sheet (FVS)), Lan (ref. 26), Lamar (ref. 27) (called Vortex-lattice Method bas- 
ed on Suction Analogy (VLM-SA)) and with data. Numerical characteristics 
and results of the present method will be presented and discussed next. 

4.1 Numerical Characteristics 

The numerical results predicted by the present method have the following 
numerical characteristics: 

(a) Convergence Characteristics 

In the present method, numerical convergence can be expected during 
the iteration process if certain restrictions are imposed on the movements 
of the free sheet and vortex core (see sections 3.5 and 3.8), 

(b) Initial Vortex Shape 

For the angle of attack below 20 degrees, the starting solution is 
taken to be the initial vortex shape from the 20 degrees angle of attack 
solution. 

(c) Spanwise Strips 

In reference 18, Mehrotra made a parametric study to find a relation 
between the aspect ratio and the number of spanwise strips which should be 
used to represent the lifting pressure on delta wings in order to obtain 
reasonably accurate aerodynamic results. The relation just mentioned is re- 
produced and displayed graphically in Fig. 10. From this figure eight 
vortex strips over the wing semi-span are seen to be needed for simple 
deltas in the aspect ratio range from 1.0 to 1.75. Because of this study 
and because of the strake-wing and double delta configurations to be treated 
here being also of low aspect ratio, the same number was used initially. 
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Due to the leading-edge kink these configurations had to be represented by 
two spanwise regions per semi-span. Experience with the QVLM of Lan’s 
(ref. 19) has shown that for two spanwise region solutions, it is better 
to have a different number of spanwise strips in each region. The final 
choice of spanwise strips used here is four for the strake region and five 
for the wing region. 

(d) Chordwise Vortex Elements 

All results to be presented have been calculated by using six chordwise 
vortex elements for all configurations. 

(e) Reference Lift Coefficient 

The reference lift coefficient discussed previously (see section 3.8) 
for use in adjusting the core movement is based on the suction analogy as 
implemented by Lan (ref. 26). 

(f) Complete Leading-Edge Separation 

In this study, only wings having complete leading-edge separation will 
be considered. 

4.2 Numerical Results 

This section presents numerical results by the present method based 
on the preceding ideas. These results are compared with other methods 

and data for certain test configurations which are shown in Figs. 11 to 15. 
Note that they consist of a flat delta wing (Fig. 11) , a conical cambered 
delta wing (Fig. 12), a delta wing with leading-edge vortex flap (Fig. 13), 
a double delta wing (Fig. 14) and a strake wing (Fig. 15). The one vortex 
core system is used for the delta wings while the two vortex-core system is 
used for the others. 


Overall Forces and Moments 


The overall aerodynamic characteristics are presented in Figs. 16 to 20. 
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All theoretical curves for the results generated are graphed by fairing 
through those at the a's given in the following table. 


Method 

Solution a's 

Present method 


Flat delta wing 

5°, 10°, 20.4°, 30.7° 

Conical cambered delta 

5°, 10°, 20.4°, 30.7° 

Delta wing with LEVF 


Double delta wing 


Strake Wing 


Mehrotra 

5°, 10°, 15°, 20°, 25°, 30° 

Lan 

All a's 

VLM-SA 

All a's 

FVS 

10°, 15°, 20°, 25°, 30° 


Delta Wing 

(a) Planar 

The results for the flat delta wing are shown in Fig. 16. It is seen 

that and can be accurately predicted by the present method as well as 
i 

by the other methods. However, the present method predicts the to be 
more positive (i.e. more nose up) at high angles of attack, whereas the 
methods of VLM-SA and Lan predict more nose down pitching moment. This 
nose up moment of the present method could be due to the estimated peak 
load being too forward. Mehrotra’s results have good agreement with data 
(ref. 22) even though he used a diffused vortex model. 

(b) Conical camber 

The results for a conical cambered delta wing are presented in Fig. 17. 
The results from ban's method have the best agreement with data. The 
results are underpredicted by the present method at lower angles of attack 
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and by the VLM-SA at higher angles of attack. FVS results (ref. 30) are 
generally low. The C results of Fig. 17(b) show that the FVS method is 
best at angles of attack below 20 degrees and the present method is best 
for angles of attack above 20 degrees. The other methods all predict 
values of C_ which are too low over the angle of attack range. Similar 
to the case for the flat delta, the present method predicts a pitch up 
tendency for the cambered wing. 

(c) Leading-edge vortex flap 

Recently, there has been much interest in applying the leading-edge 
vortex flap to increase the lift-drag ratio under maneuvering conditions. 

Some low speed results for a 74°-delta wing have been reported in reference 29. 
The geometry of the planform is shown in Fig. 13 and it was modeled theore- 
tically with its undeflected tralling-edge flap omitted for simplicity. The 
results are presented in Fig. 18 where the Cjj and have been based 
on the total planform area excluding the leading-edge flap. It is seen 
that the predicted results by the methods of Lan and VLM-SA agree well with 
data. The C and results estimated by the presented method slightly 

\a 

exceed the values of the data at moderate angles of attack. The C re- 

TQ 

suits show the present method to have a more nose up moment than the data 
for angles of attack above 12 degrees. 

Having compared the results in the one vortex-core system, it is of 
Interest to see the aerodynamic characteristics in the two vortex-core system. 

Double Delta Wing 

For the double delta wing (ref. 28), the comparison is made in Fig. 19. 
The results show that the and are well predicted by the present 

1j U • 

1 

method and the methods of Lan and VLM-SA. At high angles of attack the 
present method predicts the results better than others. 
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Strake Wing 


From reference 28, it is understood that the aerodynamic advantage of 
strake-wing planforms are (1) the transonic maneuvering capability, and 
(2) the utilization of vortex lift. For the strake-wing configuration of 
Fig. 15 (ref. 31), the results are presented in Fig. 20. It is seen that 
the predicted by the present method has good agreement with data whereas 
the method of Lan overpredicts the data in the high angle of lattack region. 
The experimental data in Fig. 20(b) shows that the pitching moments become 
more positive as the angle of attack is increased. This could be due to the 
wing vortex breakdown on the main wing. From the geometry of strake wing in 
Fig. 15, it can be seen that the wing area is much larger than the strake. 

As shown in the data in Fig. 20(a), the values start to decrease at about 

14 degrees of angle of attack, signifying the start of the wing vortex 
bursting. 


Lifting Pressure Distribution 

Figs. 21 through 29 show the spanwise pressure distributions for different 
planforms at specified angles of attack and for constant x- locations. The 
results are discussed below. 


Delta Wing 

(a) Planar 

For the flat delta wing, the spanwise pressure distributions are cal- 
culated at angles of attack equal to 10.2°, 20.4° and 30.7°, and compared 
with data at the two x-locations, x/c^^ = 0.5333 and 0.9333 in Figs. 21 to 23. 
In general, the pressure peak obtained by using Mehrotra's method is lower 
than the data near midchord, higher near the trailing-edge and is shifted 
toward the root chord. One objective in the present method is to improve 
the pressure distribution of Mehrotra's method, so as to allow for the occur- 
rence of a sharper pressure peak. The predicted pressure distributions 
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of the present method more closely resemble the data than do those of Meh- 
rotra. At the low angle of attack (a = 10.2°), the present method predicts 
a pressure peak that is inboard of the data at midchord and is lower near the 
trailing-edge. This may be due to the free sheet being too large, since its 
initial shape and core location is taken from the 20-degree solution and the 
size of free sheet will not shrink in the iteration process. Thus, the in- 
board pressure distribution are higher than those measured. 

(b) Conical camber 

The results of the estimated spanwise pressure distribution for the con- 
ical cambered delta are presented in Figs. 24 through 26, and are compared 
with data (ref. 22) and the FVS solution (ref. 30). Again, the presented 
results are at the angles of attack equal to 10.2°, 20.4° and 30.7° and 
comparisons are made at two x- locations, x/c„ = 0.67 and 0.935. At low and 
moderate angles of attack, the present method predicts values of the pressure 
peak which are lower than data, whereas the FVS solution predicts higher peaks. 
At high angles of attack the spanwise pressure distribution predicted by 
the present method has better agreement with data than does the FVS solution. 
There are two pressure peaks near the trailing-edge at all angles of attack. 

The Inboard peak is expected and is typical for flat delta configurations. 
However, the reason for the outboard peak is not conclusively known. Some 
possible reasons for its occurrence are as follows; (1) the wing vortex 
distribution is not located on the highly cambered surface; instead, it is 
on a flat surface; (2) the large spanwise camber slope term (9z/8y = -1.4402) 
near leading-edge is not included in the boundary condition; (3) eight span- 
wise strips may not be sufficient to resolve the lifting pressures in the tip 
region of this conically cambered wing. 

(c) Leading-edge vortex flaps 

Spanwise pressure distributions predicted by the present method at 24-de- 
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gree angle of attack and two x-locatlons, x/c„ = 0.5197 and 0.8937, are pre- 
seated in Fig. 27 for a 74°-delta wing with leading-edge vortex flap. Since 
there are no data or other theoretical results available, only the present 
theoretical results are presented on each figure. Fig. 27(a) shows that 

the pressure peak is at about 60 percent of the local semispan and that 
matches the vortex core location. Near the traillng-edge the interaction 
of the vortex core and leading-edge vortex flap produces large lifting 
pressure in the vortex flap region as seen in Fig. 27(b). 

Double Delta Wing 

From Fig. 28, the spanwise pressure distribution at 20-degree angle of 
attack at two x-locatlons, x/Cj^ = 0.5686 and 0.8099 are presented. In 
Fig. 28(a), the results near middle root chord Indicate that there are two 
pressure peaks present. The lower pressure peak is behind the highly 
swept inboard region with the higher pressure peak being on the outboard 
wing region. The two vortex cores induce large sidewash in the tip region 
which result in a large pressure peak there. Fig. 28(b) shows the lifting 
pressures near the trailing-edge and from this figure it is seen that another 
pressure peak has been produced. This outboard peak may be due to the large 
sidewash induced near the tip due to the interacting inboard-located vortex 
cores. The exact reasons for this require additional investigation. 

St rake Wing 

From Fig. 29, the spanwise pressure distribution at an angle of attack of 

12 degrees at three constant chord stations, x/c = 0.4,0. 7 and 0.9, are 

presented. In Fig. 29(a), the results near the kink region show that the lifting 

pressure are mainly from strake vortex. In Figs. 29(a) to 29(c), the predicted 
pressure distributions are generally higher than experimental values, the 

large pressure peak near the tip region of each figure is mainly due to the 
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effect of side-edge vortex system. Near the trailing edge the flow field is 
dominated by the wing with most of the predicted lifting pressure being on 
the wing. It should be noted that at an angle of attack of about 14° , the 
wing vortex breakdown may have occurred (ref. 37). The effect of vortex 
breakdown is not modelled in the present method. 

For a slender delta wing, a typical computer CPU time by the present 
method is 3230 seconds in 7 iterations with CDC Cyber 175 at Langley Research 
Center. This is to be compared with 4540 seconds in 10 iterations by Mehrotra's 
model. For the double delta wing analyzed in this report, the CPU time is 
5620 seconds in 8 iterations as compared with 22,450 seconds in 8 iterations 
by the Boeing code (Ref. 14). If the first two iterations in Mehrotra's model 
can be by-passed by directly assuming a reasonable initial configuration, the 
saving in CPU time is about 1200 seconds for simple delta configurations. 

4.3 Mean Square Error Computation. 

To check the accuracy of the results predicted by the present and other 
theoretical methods, the so-called mean square error (MSE) will be used. 

Mean square error is defined as : 

MSE = - d^)^)/K (111) 

where d^ is the i th experimental data value, K the total number of data points, 

and d^ the interpolated theoretical results at the location of data point d^. 

Thus, the typical results of mean square error of each theoretical method to the 

overall lift coefficient and lifting pressure distributions are as shown in the 

following: 

MSE for Overall Lift Coefficient 
(a) a flat delta wing (Fig. 16(a)) 

Present Theory 0.0209 

Lan's Method 0.0172 

VLM-SA 0.0045 


Mehrotra's Method 


0.0242 
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(b) a conically cambered wing (Fig. 17(a)) 


Present Theory 
Lan's Method 
VLM-SA 
FVS Method 

(c) a delta wing with L. 
Present Theory 
Lan's Method 
VLM-SA 

(d) a double delta wLng 
Present Theory 
Lan's Method 
VLM-SA 


0.0407 
0.0237 
0.0423 
0.0470 

. vortex flap (Fig. 18(a)) 
0.0405 
0.0378 
0.0628 
(Fig. 19 (a)) 

0.0351 

0.0330 

0.0534 


(e) a strake-wing configuration (Fig. 20 (a)) 
Present Theory 0.0724 

Lan's Method 0.1197 
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MSE for Lifting Pressure Distribution 


Configuration 


a, deg. 

Present Method 

Mehrotra 

Flat Delta, Fig. 21(a) 0.5333 

10.2 

0.1694 

0.1198 

" " Fig. 21(b) 

0.9333 

10.2 

0.1391 

0.1539 

" " Fig. 22(a) 

0.5333 

20.4 

0.2642 

0.3804 

" " Fig. 22(b) 

0.9333 

20.4 

0.1060 

0.1906 

" " Fig. 23(a) 

0.5333 

30.7 

0.5011 

0.7773 

" " Fig. 23(b) 

0.9333 

30.7 

0.1496 

0.4111 

Configuration 

x/c^ 

“» deg. 

Present Method 

FVS Method 

Delta with Conical 
Camber, Fig. 24(a) 

0.67 

10.2 

0.1540 

0.1744 

" Fig. 24(b) 

0.935 

10.2 

0.1470 

0.2163 

" Fig. 25(a) 

0.67 

20.4 

0.3146 

0.3027 

" Fig. 25(b) 

0.935 

20.4 

0.1690 

0.1358 

" Fig. 26(a) 

0.67 

30.7 

0.5090 

0.4868 

" Fig. 26(b) 

0.935 

30.7 

0.1491 

0.3139 


For the flat delta wing, the mean square errors for the lifting 
pressure indicate that the present theory is better than Mehrotra's model. 
For a delta wing with the conical camber, the present results are also good. 
However, the lower values of MSE of present results partially may be due 
to the appearance of an unexpected pressure peak near the tip region. 
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5 . Concluding Remarks 

A free sheet-vortex core model to predict pressure loadings and 
overall aerodynamic characteristics of low aspect-ratio wings with edge 
vortex separation has been developed, by incorporating a vortex core in 
the free sheet of Mehrotra's method. In the present method, Mehrotra's 
model is used in the first two iterations to establish the vortex fila- 
ment system. Based on Smith's empirical results, the initial vortex 
core location is adjusted to Improve the starting solution. Introduc- 
tion of a vortex core in the flow model increases the numerical diffi- 
culty to obtain a converged solution. The problem is solved by using 
different sets of relaxation parameters for and by imposing certain 
restrictions on movements of the free sheet and the vortex core. For 
low to moderate angles of attack ( a < 20® ) , the starting solution is 
taken to be the Initial shape of the 20-degree angle of attack solution. 

Comparison with available data indicates that: (1) the present 

method is generally accurate in predicting the lift and induced drag coeff- 
icients but the predicted pitching moment is too positive; (2) the 
spanwlse lifting pressure distributions estimated by the one vortex core 

solution of the present method are significantly better than the results 

of Mehrotra's method relative to the pressure peak values for the flat 

delta; (3) the two vortex-core system applied to the double delta and 

strake-wing produces overall aerodjmamic characteristics which have good 

agreement with data except for the pitching moment ; (4) the computer 

time for the present method is about two thirds of that of Mehrotra's 

method for a delta wing and one quarter of that of FVS method for a 

double delta wing configuration. 

In view of the Improvements embodied in the present method and 
possible further refinement, it would be of interest to: 
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(1) model the fuselage portion as a potential flow region, with 
the assumption of attached flow in the fuselage region and a 
very small thrust coefficient being assumed in that region; 

(2) examine the fuselage effect on wing vortex flow; 

(3) examine the canard effect on wing vortex flow; 

(4) compare the chordwise pressure distribution with available 
data to improve the pitching moment computation. 
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Appendix A 


An Improved Lifting Pressure Expression 


From reference 18, it is understood that the total lifting pressure, for 

the wing with leading-edge vortex element system, is contributed by wing 

bound elements, wing strearawise vortices, leading-edge bound elements and 

the strearawise vortices of the leading-edge vortex element system. Hence, 

til ttl 

the total lifting pressure at the i bound element of the j strip control 
station has the following expression: 


AC = (AC ) + (AC ) 

p . . p . , T p . . B 

J.i J.i l.i 


(A.l) 


where (AC is the lifting pressure due to streamwlse vortices; (AC )g 

^ j , i ^3,1 

the lifting pressure due to bound elements. 


A.l Derivation of (AC )„ 

From Fig. Al, it is clear that the lifting pressure due to strearawise 
vortices of wing and leading-edge vortex eleraent are calculated along the 
common edge y^ of two adjoined spanwise strips. First, the force acting on 
the chordwise element of length Ax of the leading-edge vortex system is : 


F, = p V (r. V.) Ax 
i ” 3 i 

where p is the fluid density, the free stream velocity, the vortex 

strength of the leading-edge vortex element, and v^ the sidewash at 

til 

point i. Thus, the force acting at the i point on the common edge per 


(A. 2) 


unit dynamic pressure and length is: 


(-7-) = 2 r. V. 

qAx-'i j 1 

where q is the dynamic pressure. Similarly, the force acting at the same 
point due to the right leg of the strip, per unit dynamic pressure and 


(A.3) 


length, can be written as: 

A 




ydx 


(A.4) 
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Figure Al.- The leading-edge bound vortex element and wing strearawise 
vortices at the common edge 
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where y is the wing bound vortex density, the leading-edge x-coordinate of 
the trailing leg under consideration. Using the transformation, 

+ (1 - Cos 0) 


X = x„ 


(A. 5) 


the equation (A. 4) can be reduced to the form: 




ySin 6 d0 


= -V. c. / 

i ^ ^ I 

■rrc, V. p i-1 

Sln6^J 


(A.6) 


where c^ is the local chord length of the common edge at , N the number of 
bound elements in chordwise direction and 0j^ = integral has 

been reduced to a finite sum through the conventional trapezoidal rule. 


th 

For the left leg of the (j+1) strip: 


77C. V. p i-1 


(A. 7) 
.th 


Therefore, the force per unit dynamic pressure and per unit length at the 1 
point along the common edge is: 

Thus, the lifting pressure at 1^^ point along the common edge other than 
the last one has the following expression: 

^T 

(AC )' = 2 (-^) / (y... - y.) 

J > 1 J- 


(A. 9) 


When the lifting pressure is calculated at the last common edge in the span- 
wise direction, the equation has the following form: 


(AC 


j,i 




j+1 


(A. 10) 


where yj+2’ ^j+1 y- location of common edges. Note that 

equations (A. 9) and (A. 10) are evaluated at all endpoints of wing bound 
elements, i.e. along the common edge. Surface spline interpolation (ref. 32) 
is performed to obtain the lifting pressure at the control station of each 
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vortex strip. Therefore, (AC ) , the lifting pressure due to streainwise 

^j,i T 

vortices at the i bound element of the j strip control station is obtained. 


A. 2 Derivation of (AC ) 


B 


til 

From Fig. A2, the normal force per unit length acting at the i bound 


til 

element of the j strip is given by: 




= P (u^ Yy - Y^) AS^ 

2 

p ^oo 'i'l “ Sin Il>^) Ay^/Cos ip^ 

= 2 q (Uj. - v^ Tan Yj. Ay^ 


where , 

Y^ is the bound vortex density, 

Uj^ and v^ are the x and y components of the induced velocities, 
is the sweep angle of the bound element, 

|»t_ 

Ayj is the width of the j strip, 

As^ is the length of bound element, 

Y and y are the streamwise and spanwise components of the bound vortex 
X y 

density. 

til 

Thus, the lifting pressure due to the i bound element of the j strip 


IS! 


(AC 


B, . 

) = — 7 -^ = 2(u. - V. Tan ip.) y. 
P. . T, q Ay. 1 X X 'x 


(A. 11) 


For a cambered wing without dihedral, additional lifting pressure will 

be generated from the interaction of a free stream component with the stream- 

wise vortex density y^^ (ref. 26), as shown in Fig. A3. Adding this component 

of (-2 Y Sin a Sin (j) ) to the lifting pressure due to the wing bound element, 
X y 

equation (A. 11) is modified to: 
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(A. 12) 


(AC, 


■j .i B 


) =2 Cos (j)y + Sin a Sin <|)y ) Tan 


where <f> may be defined as the local dihedral angle. Tan (3z /3y), and 

y C 

may be different everywhere on a cambered wing. The correction term in 
equation (A. 12) is applied only to the wing bound element and leading-edge 
bound element so that the local dihedral effect are Included. 

At positive a, normally the wing bound elements produce positive lift- 
ing pressure near the leading edge; while the leading-edge bound elements 
produce negative lifting pressure. Thus, the leading-edge Kutta condition 
can be satisfied when the net lifting pressure is zero at the leading- 
edge. In order to substract the lifting pressure induced by the leading- 
edge bound element from that produced by the wing bounded vortex elements, 
the following prodecures will be used: 

(1) Extrapolate lifting pressure due to the wing bound elements in 
each chordwise direction to obtain the lifting pressure at the 
location of the leading-edge bound element EF. (See Fig. 3) 

Using the Fourier series relation, the lifting pressure at any 
chordwise location is determined by the following equation: 

AC Sin 6 = a + a„ Cos £6 

p o ^ £ 


(A. 13) 


£=1 


where SlnG is included to eliminate the known square root singularity of the 
lifting pressure at the leading - and traillng-edges. 

Therefore, 


— { AC Sine de 

TT ^ P 

0 


Oi 


1 N 
N Z 
k=l 


AC Sin 9. 

Pk 


2 n 2 N 

^£=“ ^0 Cos£6 d0 ACp^SinSj^ Cos£0j^ 

0 = N 

k 2N » 


(A. 14) 


53 



where N is the number of wing bound elements in the chordwise 
direction. The integrals for' Fourier coefficients are re- 
duced to finite sums through the mid-point trapezoidal rule. 

(2) Subtract suction pressure induced by the leading-edge bound 
elements. First, the concentrated strength of the leading- 
edge vortex element is converted to the vortex density. Since, 


r = / Y dx 

(A. 15) 

from the transformation. 


c . 

X = Xjj^ + (1 - Cos 0) 

(A. 16) 

it can be obtained that 


c . 

Fj = ^ /Y Sin 0 d0 

(A. 17) 


It is assumed that the concentrated vorticity due to the leading- 
edge vortex system is distributed near the leading-edge only. Using 
the mid- point trapezoidal rule, equation (A. 17) becomes: 
i: c. Sin 6, 


r =- 1 Y 

j 2(N+1) 


(A.18) 


or. 


2(N + 1) T. 
IT c. Sin 9. 


(A. 19) 


where 0, = 


2(N + !)• 


Thus, the decrease of lifting pressure 

.th 


at the leading-edge bound element of the j strip is: 


^ “^decrease 


-2 



(v. Cos (() + Sin a Sin 4) ) Tan 

1 y y 



where the subscript £e means that u., v., 6 and ib . 

J j’ y 

at the leading-edge bound vortex element of the 
the actual lifting pressure at the location of the 
bound element of the strip is: 


e 

are evaluated 
strip. Hence, 
leading-edge 


(A. 20) 
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(1C ) = fdc^ )- + (1C„ )-] + (1C„), 

"j.l '■ ‘’J.i® "3,1 lie 


(A.21) 


decrease 


The lifting pressures defined by equation (A.l) at a discrete number of 
points are Fourier-analyzed through equation (A. 14) so that lifting pressures 
at any locations can be calculated. 
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Appendix B 


Correction to Sidewash Calculation in Discrete Vortex Method 


In developing the free vortex filament theory for calculating aerodynamic 
characteristics of wings with edge-separated vortex flow, the position of the 
concentrated vortex core is adjusted in the direction of decreasing Kutta- 
Joukowsky force. The magnitude of the force, and hence the new position of the 
vortex core, depends on the calculated induced velocity components on the core. 
It is well known that the induced sidewash cannot be accurately calculated by 
the discrete vortex approximation at points close to the vortex sheet. There- 
fore, the vortex core location and the wing aerodynamic characteristics will 
not be accurately evaluated. A method to improve the accuracy of calculating 
sidewash by discrete vortex approximation will be detailed in the following 
paragraphs. 


B.l A Test Case 


Assume a vortex distribution with constant density (to = lo T) exist in a 
rectangular region, as shown in Fig. Bl. The induced velocity vector at any 
point in space with position vector is then given by: 


g. ^ - R) X (0 


= TT / 


2 

4t 7 ■' s 


d s 




= (4 - z)j - (n - y)k 

[(C-x)2 + b2 (n-y)" + B^ (4-z)2]3/2 


d5 dn (B.l) 


where 



2 2 2 2 2 2 
Rg = (5 - x)^ + 6 (n - y) + B (? - z)^ 

Rj^ - ^ = (5 - x)? + (n - y)l + (C - z)^ 

Since the vortex density to is assumed to be constant over the entire region, 

equation (B.l) can be integrated (ref. 33) to give a closed-form expression 

for the sidewash V as: 
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Figure Bl.- Vortex distribution in a rectangular region 


N,n 



\ / 


\< 

Q(x,y,z) 


21 and 34 are vortex 
segments of a strip 


Figure B2.- Local panel geometry 
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(y - y) (a - x) 


„ 0)3 )„ -1 
v=^^Tan 


- Tan 


? - z 


[(a - x)^ +3^ (y - y)^ + 3^ (C - z)^ J 


2-]h 


+ Tan 


- Tan 


■1 

(-y)(a - x) 



1 

N 

[(a - x)^ + 3^ y^ + 3^ (4 - ; 
(| - y) (x) 


H 

U-^l| 

1 

^x^ + 3^ (| - y)^ + 3^ (C - 

(-X) (y) 

z)^ 


U - z| 1 

[x^ + 3^ y^ + 3^ (4 - z)^ 

h 

(B.2) 


Equation (B.2) will be used to check the accuracy of approximate methods of cal- 
culating the Integrals in equation (B.l). 


B.2 Approximate Methods 

The chordwise integration in equation (B.l) is performed exactly so that 
the sidewash component becomes: 


4 IT 


o)(n) d n 


a - X 


o (q - y) + (C - z) [ jT(a - x)^ + 3^ (q - y)^ 

(B.3) 


+ 3^ (C - z)^ p + 3^ (q - y)^ + 3^ (C - z)^ ]^J 

In conventional discretized approximation, equation (B.3) is reduced to: 




“k 


a - X 


k=l (tIj^ - y)^ + (? - z) [ |^(a - x) +3 (qj^ - y)‘ 


X 


+ 3^ (4 - zfp [x^ + 6^ (qj^ - yf + 3^ (C - 


(B.4) 


The results in equation (B.4) are equivalent to applying mid-point trapezoidal 
rule to evaluating the spanwise integral in equation (B.3). However, when 
(^ - z) is small, the integrand in equation (B.3) contains a second-order sing- 
ularity at q = y, so that equation (B.4) is expected to be inaccurate. To 
improve equation (B.4) and yet keep the method as simple as possible, the 
following consideration is applied. 
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For convenience, let 


f(n) = 


a - X 


jj(a - x)^ + 6^ (n - y)^ + 3^ (? - z)^ 


X 


[x^ + 3^ (n - y)^ + 3^ (? - z)^ 

If < y < equation (B.3) can be rewritten as: 


^ (g - z)3' 
4x 


f( /’^i+Z 

L o 


\+l b/2 

+ / ) 


^ (n) f (n) 


li+i (n - y) 


j dJ 

+ (? - Z )^ J 


The first and the last integrals can be reduced to finite sums as usual: 

g - ^ , , 

Li , 1 / n 2 

^k=l (tIj^ - y) + (C - z) 


where 


V = 


“i ^ ^^i^ ^’^i 


(B.5) 


(B.6) 


(n^ - y)^ + (c - z)^ 


+ % 


“i+1 ^ 


i+1^ “"i+1 


N 

+ z 


“k ^ 


(rii^l “ + (? - z)^ k=i+l ir\^ - y)^ + 


» + e1 (B.7) 

(C - z)^ J 


1+1 

E = / -^ oi (n) f (n) 

^^i (n - y)^ + (C - z)^ 


dn 


= Lti) -- y( 2 ) f . (y) dn + 0 )(y) f (y) ^ 


(n - y)^ + (? - z)^ 


’^i (n - y)^ + (c - z)' 


u(n^) f (n^) - U) (y) f (y)^ “(’^i+1^ ^ ^’^i+1^ - o) (y) f (y) 

= u ' + ^ 


— *2 


(n^ - y)^ + (4 - z)^ 


-An 


^ c^xIXCzl fran-1 - Tan-^ 

l?-zl L 4-z 


(hi +1 - y) + (c - z)‘ 
n. - y 


i+l 


(B.8) 


In equation (B.8), the first integral was evaluated by conventional trapezoidal 
rule which is probably not accurate in this case. However, this is necessary 
to keep the method simple so that it is easily extendable to more general 
applications. A correction factor may be applied, as will be indicated later 
(see equation (B.IO)). 


59 



Substituting equation (B.8) into equation (B.7) results in: 


V = 


- z)S 


2p N oiy, f (Hj^) An. 

^ 2 

^k=l (n,_ - y) + 


Att 


- 


u(y) f (y) An. 


(\ - y) + (c - z) 
w (y) f (y) An^^j^ 


(n. - y) + (? - z)' 


(Hi+l - y) + (c - z) 


1 , £0 (y) f (y) 3 


I^Tan"^ 


\+i - y 

? - z 


- Tan 


-1 \ - y ] 

z| J 


(B.9) 


The first term in equation (B.9) represents exactly what would be given by 
the conventional discretized approximation (see equation (B.4)). Therefore, 
the remaining terms represent a correction C: 

t»(y) f(y) An. 


G ,(c - z)sM“W N 

“ 2 '■ 4 t: M , .2 ^ 

L - y) + (C - z) 


■i+1 


(Hi^l - y)^ + (C - z)^J 


^2 I , n,.j_T - y 


+ ^ ) <Tan-l -4±i r - Tan 

?-z ' 47T ■^j C-Z 


•1 \ - y 

? - Z I 


(B.IO) 


where G is a correction factor mentioned above and is taken to be 1.1 in the 
following numerical study. 


B.3 General Applications 

To apply equation (B.IO) to general situations involving such as leading- 
edge vortex sheet, trailing-vortex mutual Interaction, etc., all geometric 
quantities involved in equation (B.IO) should be evaluated based on local 
panel coordinate system. To illustrate this, consider two vortex segments 
21 and 34 in Fig. B2 and "P" is a point on a plane normal to two diagonal 
vectors (T^^ and T 2 ) and is obtained by projecting the control point Q along 


the normal vector. The vectors Tj^ and T 2 can be calculated as: 
^1 ^^3 “ ^ 1 ) ^ + (y 3 - y^) 3 + (Z 3 - Zj^) ic 

^2 " ^^4 - ^2^ ^ ^y4 " ^2^ ^ ^^4 “ ^2^ ^ 


(B.ll) 

(B.12) 
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It follows that a normal vector N may be defined as; 


N = X T2 = 


-i- 

1 


X3 Xj_ 


- X2 


->• 

3 

73 - 


t 


Z3 - 


^4 - ^2 

= i [(73 - y^) (Z4 - 2^) - (Z3 - z^) (y^ - y2)j 

+ j |^(Z3 - z^) (x^ - X2) - (X3 - Xj^) (z^ - Z2>j 

+ t [(X3 - Xj^) (y^ - y^) - (y^ - y^) (x^ - X2)J 


Z4 - ^2 


= N i + N j+N k 
X y z 


(B.13) 


A unit vector associated with N can be defined as: 
N 


n = 


1^1 


(B.14) 


To determine the coordinate of "P", the distance PQ is needed. However, 

PQ is just the projection of vector 2Q onto the normal vector n : 

PQ = 2^ • n 

= (x - X2) ~ ^2^ "y “ ^2^ "z (B.15) 

It follows that the coordinates of point P are: 

Xp = (x - ^2^"x 
yp = (y - y2)ny 


Zp = (z - Z2)n^ 


(B.16) 

To find the distances "x" and "a - x" in equation (B.IO) (through f(y)). 


the unit vectors and 03^ are needed: 


"21 


‘34 


2l% 


(xj^ - X 2 ) 1 + (yj^ - y2> 3 + (z^ - ^2^ ^ 

[<Xi - + (yi - y2>^ + ^^1 - ^2>X 

(x^ - X3) t + (y^ - y3> ^ + (z^ - Z3) k 


[(x^ - X3) + (y^ - y3> + (z^ - Z3) ]' 

Then, the "x" distance in equation (B.IO) can be replaced by: 


(B.17) 


(B.18) 


^ ^^2p‘ “21 ^3p* “34^ 


(B.19) 
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where 


■§2p = (Xp - X2) 1 + (Yp - y2> J + (Zp - Z2> ^ (B.20) 

S^p = (Xp - x^) i + (Yp - 73 ) j + (Zp - z^) k (B.21) 

On the other hand, the "a - x" distance In equation (B.IO) can be replaced by; 

4 " ^ ^^Ip* ^21 ^34^ 

where 

? 3 p = (x^ - Xp) 1 + (y^ - Yp) 1 + (Z 3 - Zp) t (B.23) 

^4p = (^ 4 " ^p> ^ + ^^4 " yp> J + (^4 " ^p^ ^ 

To determine "y - and - y" in equation (B.IO), cross products 

of S 2 p with and S^p with u^^ are needed. Hence, "y - will be replaced 

by: 

= 1 ^ 2 p ^ ^21 1 

Similarly, - y" will be replaced by; 

Yjj^ “ ^^3p ^ ^^34^ (B.26) 

Finally, "? - z" in equation (B.IO) must be determined. It is simply 

replaced by -PQ in equation (B.16): 

^ - y = - [(Xp - x)^ + (Yp - y)^ + (Zp - z)^ (B.27) 

B.4 Numerical Results 

The locations of control points chosen for illustration are indicated 
in Fig. B3. The results are compared in Figs. B4 through B7 for M=0. Except 
for the case where the control point is directly under one, of the discrete 
vortices (i.e. point B) , all results with the derived correction terms "C" 
in equation (B.IO) appear to be good. In applications to free vortex-filament 
program, the control points on the concentrated vortex core are chosen to be 
nearly at the middle of vortex strips formed by the leading-edge free vortex 
filaments. 
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Figure B5 








Appendix C 

Determining Initial Core Location by Vector Analysis 

After exercising Mehrotra's code (ref. 18) through two iterations and 
establishing the entire leading-edge vortex element system, the centroid of 
the leading-edge vortex filament system is determined from equations (81) and 
(82). Therefore, the initial configuration of leading-edge vortex filament 
system and its centroid are as shown in Fig. Cl. 

Based on the vector analysis of the cross product, the procedures of 
finding the intersection between the leading-edge vortex filament system and 
its line of centroid are discussed in the following. 

Note that each leading-edge vortex filament consists of a series of 
short straight segments. It is desired to find the intersection point between 
the leading-edge vortex element system and its line of centroid on the X-Y plane. 
Thus, only x and y coordinates of each segment will be considered in the cross 
product analysis. Assume that ^ segment along the leading-edge vortex 

filament system and is a segment along the line of centroid. Then, 
and have an intersection if the following conditions are satisfied: 

(1) Resulting vectors of x ^4 and 8^2 ^ are in the opposite 
direction. 

-y “V -V “>■ 

(2) Resulting vectors of ^ ^14 ^ ^12 opposite 

direction. 

For example, in Fig. C2, by using the right hand rule and cross product analysis, 
^^2 ^34 have an intersection point. Similarly, ^^^23^^^ ^34 have no inter- 

sections in Fig. C3. 
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X 


Figure Cl.- Leading-edge vortex filament system and its line of 
centroid 
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Appendix D 


Corrections to Induced Velocities Due to Vortex Core 


It is clear from Appendix A that the lifting pressure at the wing bound 
element will not be predicted accurately, if the computation of sidewash due 
to the vortex core is not treated properly. In fact, a very large sidewash 
would result whenever the vortex core is too close to the wing surface. To 
Improve the sidewash computation due to the vortex core, and hence, the 
lifting pressure and the aerodynamic characteristics, the computed sidewash 
due to the vortex core will be corrected as follows. 

In reference 34 (on page 592) , Sir Horace Lamb showed that the circum- 
ferential velocity around a circle with radius r due to a line vortex of 
strength F is: 


"^0 2irr 


(1 - 


(D.l) 


Later, Kutler in reference 35 defined the approximate radius as: 


r = 4vt- 
c o 


(D.2) 


Thus, the tangential velocity in equation (D.l) becomes: 


V 


0 


r 

2irr 



2 

c) 


(D.3) 


In equation (D.l) through (D.3), r is the radial distance from the center of 
the concentrated line vortex, F the circulation strength, r^ the approximate 
core radius, v the kinematic viscosity, and t^ the initial time. 

From equation (D.3), after expanding the exponential term, the tangential 
velocity becomes: 


V 


0 


F 

2irr 



2 

c) 


F 

2nr 



(D.4) 
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The above equation is obtained from the assumption that the ratio r/r is 

c 

r 

less than one. In the present analysis, the term of is to be regarded as 
the computed sidewash without correction, and all terms in equation (D.4) are 
defined as follows. 


In Fig. Dl, the vortex core segment is on the right wing and ? 


’78 


on 


the left wing. a leading-edge vector. Point 5 is the wing bound 

element location where the induced velocities due to vortex core segments 
and are to be calculated. For the right core segment the 


78 


34’ 


characteristic core radius r^ is defined as 3/4 of the perpendicular distance 
from point 4 to line ^22’ 


3 It 12 

r = 7* S-, X 

c 4 ' 14 

■'l2l 


1 ^. 


(D.5) 


The numerical factor, 3/4, in equation (D.5) is an assumed empirical value. 
The radial distance r of the wing bound element at point 5 to the vortex 
core segment is calculated as : 


= I? 


54 ^ 


34 


^34 I 


(D.6) 


Thus, the radial distance r is the perpendicular distance from point 5 to 

Based on the assumption in equation (D.4), the computed induced side- 

wash v' is to be modified only when the ratio r/r is less than one. The 

c 

corrected Induced sidewash is: 


V = V* (— ) , r/r < 1.0 

r c — 

c 


(D.7) 


that is the location at which the Induced velocity is to be computed is inside 
the vortex core region. 

Similarly, at the right wing bound element at point 5, the sidewash due 

to the left core segment is calculated, based on the following r and r: 

/o c 
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X 


X 


(D.8) 


r 

c 






r 



78 


= I t 
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(D.9) 


Again, the correction of the computed sidewash will be applied only to those po- 
ints inside the vortex core region, and the modified sidewash has the same 
expression as that of equation (D.7). Note that the right wing vortex core 
segment and the left wing vortex core segment S^g are geometrically 
similar. Therefore, in the actual computation, r^ and r due to ^^g are 
based on the right wing vortex core segment with the bound element loca- 
tion switched from point 5 on the right wing to its image point 6 on the left 
wing. 
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Figure 1.- Simple delta wing geometry without leading-edge vortex system 
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Figure 2.- Strake wing geometry without leading-edge and side-edge 


vortex system 
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Figure 5.- Vortex segment geometry 




) Top view 


of the initial free sheet and vortex core 





(c) Side view 
Figure 6.- Concluded. 
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Figure 7.- Lifting pressure components on the wing camber surface 
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Figure 8.- Definition of n^ and t^ 



Figure 9.- Sketch of diffused vortex filament and vortex core 
models 


82 


0 


2 


10 


12 


1 * 6 8 
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Figure 10.- Variation of number of spanwise strips with 
aspect ratio for simple delta wings 
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Figure 13.- Geometry of a 74°-Delta Wing with Leading-edge 
Vortex Flap 
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Figure 14.- Geometry of a Double Delta Wing 
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Figure 20.- Longitudinal Aerodynamic Characteristics of a Strake Wing; M = 0. 
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(b) X/C = 0.8937 
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Figure 27.- Concluded. 
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(b) x/c = 0.7 
Figure 29. ~ Continued. 
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(c) x/c = 0.9 
Figure 29.- Concluded. 
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